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Abstract 

In this paper we propose a simple mean-field "toy" model for the liquid-glass phase 
transition. This is the system of N point-like particles confined in a finite volume of a D- 
dimensional space interacting via infinite-range oscillating potential. In the framework of the 
replica approach it is shown that such a system exhibits the phase transition between the 
high-temperature liquid phase and the low-temperature glass phase. This phase transition is 
described in terms of the standard one-step replica symmetry breaking scheme. 
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1 Introduction 



The problem of the liquid-glass phase transitions attracts permanent interest during last decades 
(for reviews see e.g. 0]). In the recent years in additional to the traditional experimental and 
phenomenological investigations a notable progress has been achieved in a first principle statistical 
mechanical study of the glass phase (see [[| and references therein). Leaving apart a wide scope 
of non-equilibrium properties of glasses, in a pure statistical mechanical approach one is aiming to 
investigate the thermodynamical properties of the iV-particle system with the two body interparticle 
interactions described by the Hamiltonian 



N 

H[m]= £ Ufa-,*), (i.i) 

where x, is the vector in a D-dimensional space which points the position of the i-th particle, and 
t/(x) is the interparticle potential (one, of course, can consider more general systems with two 
or more sorts of particles and different potentials for different particles). In a realistic systems the 
potential U (x) must be attractive (and sufficiently quickly decaying) at large distances and strongly 
repulsive at short distances. 

In a very simplified form the scheme of calculations is supposed to look as follows. Instead of 
the plain (and hopeless) integration over all positions of the particles in the partition function 



N 

n dD 

.1=1 



x, 



exp(-/3iy[x i ]) (1.2) 



in the systematic approach, first of all, one should find (or rather guess) the space structure and the 
energy of the (zero-temperature) ground state configuration of the Hamiltonian, and then, using 
one or another approximation, integrate over the fluctuations around this state. Depending on the 
value of the temperature one eventually finds that either the proposed non-trivial ground state 
is stable with respect to the fluctuations (at low enough temperatures), or (at sufficiently high 
temperatures) the fluctuations destroy proposed ground state and the thermodynamic state of the 
system is a liquid where all the particles are delocalized. 

Indeed, the above scheme works rather well if the low temperature ground state of the system is 
ordered. In this case we are dealing with the crystal which is characterized by one or another global 
symmetry breaking, and which can be sufficiently easily described analytically. However, in the 
study of the disordered glass state, the situation becomes much more complicated, because in this 
case the low-temperature solid state is characterized by the random positions of the particles, and 
it is impossible to define what kind of the global symmetry is broken. All that resembles a lot the 
problem we are facing in spin-glasses, where the spins are getting frozen in a random state which 
can not be characterized by any apparent global symmetry breaking. The only but quite essential 
difference is that, unlike spin-glasses, here we do not have quenched disorder installed in the initial 
Hamiltonian. Nevertheless, the ideas borrowed from the spin-glass theory, and in particular the use 
of the replica technique (which, as we realize now, has much deeper meaning than just a technical 
trick to pass over the averaging of the logarithm of the partition function) turned out to be quite 
fruitful also for structural glasses, as was shown in the series of papers by Mezard and Parisi E3J. 
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To demonstrate the effect of a symmetry breaking in ordered magnetic systems one can introduce 
a conjugated field coupled to the order parameter, which at the end (after taking the thermodynamic 
limit) is set to zero. In spin-glasses the same effect can be achieved by introducing several weakly 
coupled copies (replicas) of the original system. In a similar way, to demonstrate the effect of the 
freezing into a random glass state in the system of particles described by the Hamiltonian, eq. fll.ip , 
let us introduce two identical copies of the same system described the following Hamiltonian 

N N N 

# = E u (*i - x i) + E ^(y* - y*) + e E - » ( L3 ) 

which contains a weak (controlled by the parameter e) attractive potential W between particles x 
and y of the two systems. After taking the thermodynamic limit the parameter e must be set to 
zero. This trick can result in two types of situations: 

(1) After taking the limit e — > the particles of the two systems becomes completely independent 
(uncorrelated). This would indicate that the particles do not have a "memory" of their spatial 
positions so that they are free to move (non localized in space) and the original system is in the 
(high-temperature) liquid phase. 

(2) Even after taking the limit e — > the positions of the particles of the two systems remain 
correlated. This would indicate that the particles become localized in space so that the original 
system is in the (low-temperature) solid state. The order parameter describing this phase can be 
defined e.g. in terms of the correlation function between particles of the two copies of the system. 

In fact, like in spin glasses, to obtain more detailed information about this type of phase tran- 
sition, instead of the two copies it is more convenient to introduce a general n replicas of the 
original system (see Section 2). It should be also noted that in real calculations the introduction of 
the supplementary attractive potential between replicas is actually not necessary. It is well known 
that e.g. in the case of the paramagnetic-ferromagnetic phase transition, instead of introducing a 
conjugated field, it is sufficient just to suppose the possibility of the global symmetry breaking to 
prove its existence afterwords. In a similar way, here we are also going to admit the possibility of 
the effective space correlations among particles of different (originally non-coupled) replicas, while 
the validity of this assumption can be checked aposteriori. 

Although the above idea is very simple, actual calculations for realistic models turn out to be 
quite sophisticated (see e.g. 0). What is lacking in this field is the analog of the mean-field models 
of spin-glasses, like the SK-model |4j] or REM ||, which despite of their non-physical nature, turned 
out to be very powerful tool for understanding the nature of the spin-glass state. The aim of this 
paper is to present a very simple mean-field model of a glass, which, despite of its completely 
non-physical structure, exhibits non-trivial liquid-glass phase transition (Section 3). Moreover, the 
nature of this phase transition turns out to be of the so-called one-step replica symmetry breaking 
type (the one, which, in particular, takes place in REM), which is characterized by a finite jump of 
the order parameter (like at the first-order phase transitions) and continuous free energy function 
(as it should be at the second-order phase transition). 
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2 Replica Calculations 



The general scheme of the replica calculations (which we will follow in the next section) can be 
described as follows. First, we take n non-coupled copies of the original system. The corresponding 
replica partition function is 



Z 



n 



N n . 

n n / d D < 

i=l a=l 



exp 1-{3Y,H[<\) , (2-1) 



a=l 



where the parameter n should be set to n — 1 et the end. Next, following the standard scheme of the 
one-step replica symmetry breaking we divide n replicas into n/m groups each consisting of m 
replicas. We will suppose that particles belonging to different groups of replicas are non-correlated, 
while those of the same group could become effectively coupled. In this case the partition function, 
eq. reduces to 

71 

Z n = [Z m p , (2.2) 
and the density of the free energy of the system is then 

1 n 

f = = -^N^ (2 ' 3) 
After taking the limit n — > 1 we get: 

f(m,P) = -^\n(Z m ) (2.4) 

where m becomes now the continuous variational parameter bounded by the condition m < 1 (this 
constraint comes from the original bound m < n in the limit n — > 1). 

The assumption that the positions x" of particles in different replicas are correlated can be 
explicitly represented as follows: 

X? = X; + u? , (2.5) 

where Xj play the role of the center of mass of the replica " molecules" , and u" are the deviations 
of the particles from the center of mass. Of course, these deviations are bounded by the condition 

m 

E < = (2-6) 

o=l 

In this way the replica partition function takes the following form: 



N 



.i=l 



JO, 



N m 



. Li=l a=l 



n n / dD * a 



N 



Hm D 5 E< 



.i=i 



va=l 



exp + 



[2.7) 



a=l 



We see that now the problem becomes similar to that of statistical systems with quenched disorder 
in replica representation: according to eq.( ^.7| ), Xj's play the role of disorder parameters, while u" 
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are the dynamical variables. As usual, after averaging over the disorder parameters {xj} we will 
get the partition function Z m represented in terms of a new replica Hamiltonian if m [uf]: 



N 



nnK u 



=1 a=l 



exp(-/3tf m [<D 



(2.{ 



where the variables uf with different replicas could become coupled. 

In the final step of calculations will be looking for the extremum of the free energy, eq.( 
with respect to the continuous parameter m, bounded by the condition m < 1. As usual in the 
replica calculations, after the parameter m (which originally was introduced to describe the number 
of replicas) is taken to be smaller than one, the minimum of the free energy becomes a maximum. 
In the present case due to the constraint, eq.( |2.6| ), the number of independent replica variables u" 
is equal to (m — 1). After taking the parameter m to be smaller than one this number formally 
becomes negative which effectively changes the minimum of the free energy into the maximum (this 
situation is quite similar to that of the standard replica calculations in disordered systems ||) 

The maximum of the free energy, eq.([2.4|), is defined by the saddle-point equation 



df(m,P)/dm = (2.9) 

If the solution of this equation, m*(/3), appears to be smaller than one, then for the physical free 
energy one finds: f(/3) = /(m*(/3), (3). This situation corresponds to the solid glassy phase, and it 
will be shown to take place only at temperatures smaller than a certain critical temperature T c . 
At T > T c in the high-temperature liquid phase the formal solution of eq.(|2.9|) yields m*(/3) > 1 
(which is not allowed by the constraint m < 1) and it can be proved that at the interval < m < 1 
the maximum of the free energy f(m, (3) is achieved at m = 1. Thus, the above scenario looks quite 
similar to that of the standard one-step replica symmetry breaking in the random energy model of 
spin-glasses [|]. 



3 The model 

The model we consider in this paper is defined by the following Hamiltonian: 

i N -I N 

tf = -^I>(l*.-*,l) + ^I>? (3-D 

V 1 * J<J 1=1 

where {xj} are the coordinates of the point-like particles in the D-dimensional space and 

E/(|x|) = V2cos(|x|) (3.2) 

is the infinite-range oscillatory interparticle interaction potential. The second (Gaussian) part of 
the Hamiltonian, eq. (|3.1| ), bounds the particles to be effectively confined in the finite volume with 
the linear size L. The parameter L 3> 1 is supposed to be large compared to the oscillation period 
of the potential, eq.( |3.2|) , and it will remain finite in the thermodynamic limit N — > oo. 
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The interaction potential, eg. ( |3.2|) , yields the alternating concentric bands of positive and neg- 
ative energy. In dimensions greater than D = 1 in a system of N particles with the Hamiltonian, 
eq.( p.l|) , it creates highly complicates "interference" pattern, so that such system appears to be 
strongly "frustrated" in a sense that there exists no optimal position of the particles which would 
satisfy all interparticle interactions. Therefore, a generic solid state (if it exists) in this model must 
be a glass and not the crystal. 



The corresponding replica partition function is 



N m cl D x. a 







m N 



ex p H= EE u(K-*]\ 



a=l i<j 



m N 

££<*■ 

a=l i=l 



2L 



(3.3) 



(for convenience we have introduced here the "natural" normalization factor (2nL)~ D / 2 for the 
integration over {x^}). After the splitting of the degrees of freedom {x"} into the center of mass 
{xj} and the deviations {u"} (eqs. (|2.5|) - (|2.6|) ), we get: 



(2ttL)" 



x exp 



N m 



n n / d°< 



i=l a=l 
m N 



N 



J] m D 5 u i 



\a=l 



N 



i=l 



n h 



X; 



X 



a m JS o N n m N 

E E u{\*i - x, + <,|) - ^ e x, 2 - £ E E« N2 



a=l i<j 



i=l 



2L 



(3.4) 



a=l i=l 



where u?- = — u". Keeping only extensive in N terms (the only relevant ones in the thermody- 
namic limit iV — > oo) for the (Gaussian) averaging over {xj} we obtain: 



-N 

Y[ / d D x;exp 
.i=i 



f3m 
' ~2L 



x 2 



exp 



m N 



v EE^(lx,-x J + <.|) 

iV a= i iK j 



(2tiL\ 



ND 
2 







m N 



2 m N N 



2N 



E ££^(|x i -x, + u«.|)[/(|x fc -x i + <|) 



a,b=l i<j k<l 



(3.5) 



where 



2nL 

(5m 



■N 

Y[ / rf D Xjexp 

i=l 



0m_^2 



2L 



(3.6) 



According to the above definition, eq.( |3.6|) , one can easily prove (under condition L 1): 



tr(|x|) = o 



(3.7) 
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f/dxilMlxjl)^, 



(3.8) 



and 



*7(| Xi |)*7(| Xj + u|)~ J (|u|)5 - 
where Jq(u) is the Bessel function. In particular, for a < 1 we have: 



(3.9) 



J {u < 1) ~ 1 - -u 2 



(3.10) 



Substituting eq. (|3.5| ) into eq. (|3.4j) and using eqs.(p.7|) and ( |3.9|) for the partition function, 
eq.(|3.4|) we obtain: 



ND ND JVg(m-l) 

m 2 p 2 (2nL) 2 



( E E«) 2 + ^EE Jb(k a - ^ + u ; - % u 



X 



2L 



a=l i=l 



N m 

n n / 

|_i=l a=l J 

2 m V 

27V 



'AT / m 

IPS*? 

i=l \a=l 



X 



(3.11) 



a, 6=1 i<j 



Assuming that all the deviations {u^} are small (only under this condition the separation of the 
degrees of freedom defined in eqs.(|2~T5|)-( |2.6| ) makes physical sense), using eq.( |3.10| ), we get 



ND ND , ND(m-l) 

m 2 fj 2 (27rL) 2 



N m r 

n n / rfD < 

.i=l a=l 



N / m 

IF E< 

i=l \a=l 



X 



exp (l/3WiV - i(/3 2 m + |) £ XK) 2 ) 



Z m = exp (-(3 2 m 2 N - -NDm\n/3 - -ND(m - l)ln(l + /3mL)\ 
V 4 2 2 / 

For the free energy density, eq. (|2.4| ), we eventually obtain: 



(3.12) 



Simple Gaussian integration over {u"} (taking into account the constraints Y^ZLi u ? = 0) yields: 



(3.13) 



(3.14) 



One can easily prove that this function of the parameter m has a unique maximum at m = m, 
defined by the saddle-point equation df /dm = 0: 



At L ^> 1 the approximate solution of the above equation can be obtained explicitly: 



(3.15) 
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m*(£) ~ -V2D\nL (3.16) 

Since the function f(j3,m) is defined only at the interval < m < 1, in the case m* > 1 the 
maximum of / is achieved at m = 1. Thus, the maximum of the function f(j3,m) (0 < m < 1) 
takes place at 

^V2Dh^L = | ;for/5 >/5 c 
m*(/3) = <! (3.17) 



1 ; for /? < & 



where 



P c ~ V2D\nL (3.18) 

Correspondingly, for the physical free energy of the system, /(/?) = f(j3, m*), we obtain the following 
result: 

r §ln/3 -fe-^tf- p c )\n(l + P C L) ; for /3 > /3 C 

/OS) = (3.19) 
[§ln/5-i/3 ;for/5</5 c 

One can easily verify that direct (without replicas) calculation of the free energy 



HP) = -^ln 



N 



i=l 



exp{-/3if[xi]} (3.20) 



yields the same result as in eq. (|3.19| ) in the high-temperature region, P < P c . Thus, plain integration 



over all particles positions (as it is done in the direct calculations, eq. (|3.20|) ) correctly describes only 
the high-temperature (liquid) phase, but it completely misses the "condensation" of the particles 
into the solid (glassy) phase at P > P c . 

The order parameter which describes the correlations of the particles inside the replica " molecules" 
in the glassy phase can be defined in the standard way: 

QHK-u 6 ) 2 ) (3.21) 
where a ^h. Since all the replicas are equivalent we can also define Q as follows: 

i m 

Q = —, rrE(( ua - u6 ) 2 ) (3-22) 

Using the constraint, eq. (|2.6|) , we get: 



Q = 7 tt(K) 2 ) (3-23) 

[m — 1) 

where the replica number a and the space vector component a are arbitrary. Note that the above 
definition of the order parameter is valid only for the glassy phase at m < 1. In the liquid phase 
(m = 1) the positions of the particles are not correlated by definition. 
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Proceeding similarly to the calculation of the replica partition function, (eq. (|3.12|) ), one finds: 



<«r> 



Straightforward calculations yield: 



'Ua=lJd D U a _ 


«) 2 6 (EZi u a ) exp 


^-|(/3 2 m+|)E^i(u° 


) 2 ) 


[nr=i/^u«]5(Er=iU a )exp 









(3.24) 



(«) ) 



m — 1 



(3m((3m — L v 
Correspondingly for the order parameter, eq. ( |3.23| ), we get: 



(3.25) 



Q 



2D 



(3.26) 



Pipm-L- 1 ) 

Substituting here the saddle-point value of m = m*(/3) obtained above, eq.( |3.16|) , at L > 1 we 
eventually find: 



Q 



2D 



2D 



(3.27) 



(3V2D\nL PPc 

We see that in the whole low-temperature region at (3 > (3 C ~ \J2D\yvL ^> 1 the value of the order 
parameter Q (as well as the value of the typical deviation ((m„) 2 )) remains small, which justify the 
approximation made in calculation of the replica partition function, eq. (|3.12|) , as well as the whole 
original idea of splitting the degrees of freedom into the centers of masses of the replica "molecules" 
and the small deviations u a , eqs. (|2.5| )-(|2.6|). 



4 Discussion 



According to eq. fl3.21p the order parameter Q describes the typical correlations of the particles 
belonging to the same replica molecule. On the other hand, according to the procedure described in 
Section 2, the positions of the particles belonging to different molecules remain non-correlated, and 
the value of the the fraction of the particles bounded in molecules is controlled by the parameter 
m*(/3). According to the standard physical interpretation of the replica parameters (see e.g. 0,0), 
the fraction of particles localized in space (a kind of a "glass condensate") must be equal to (1 — 
m*(/3)) = 1 — /3 c //3, and the fraction of particles which remain free to move must be equal to 
m*([3) = p c / p. Physically the order parameter Q describes the typical value of the space fluctuations 
of the localized particles around their equilibrium positions. 

It is interesting to note that the value of the order parameter Q computed in the previous 
Section remains finite in the whole low-temperature (glassy) phase including the phase transition 
point. In particular, according to eq. (|3.27| ) at j3 = P c we have: 



Q.* m (4.1) 
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and Q(T — > 0) — > 0. Note also that in the high-temperature phase this order parameter is not 
defined: in terms of replicas it can be defined only for m 7^ 1 while at (5 < (3 C) m*(/3) = 1, and in 
physical terms it describes the fluctuation of the localized particles which are just absent at (3 < (3 C . 

Thus, here we face the typical scenario of the phase transition with one-step replica symmetry 
breaking well known after the random energy model of spin glasses f|. On one hand, in terms of 
the free energy and others thermodynamical functions this is the phase transition of the second 
order (the free energy function is continuous at T c and it has a singularity only in the second 
derivative over the temperature, see eq.( 3.19| )), while on the other hand, the order parameter which 



characterizes the low-temperature phase has a finite value right at T c (which corresponds to the 
first-order phase transitions). 

In physical terms the present phase transition can be described as follows. Above T c all the 
particles of the system are delocalized and the system is in the liquid phase. Just below T c the 
glassy phase is characterized by a small fraction (equal to 1 — T/T c ) of particles which become 
localized at random positions in space such that their thermal fluctuations around these equilibrium 
positions described by the order parameter Q have a finite (small) value right at T c , eq.( |4.1| ). At 
lowering the temperature the fraction of localized particles increases as (1 — T/T c ) while the value 
of their thermal fluctuations around the localized positions decrease as Q(T) ~ 2DT C T, eq.( |3.27|) . 
Eventually, at T — > all the particles of the system become localized in a random positions (without 
thermal fluctuations) such that the system turns into perfectly frozen glass state. 

Of course, from the physical point of view the present model has two very important defects 
which makes it not more than a "toy model" of the glass transition. First, it has completely non- 
physical long-range structure of the interactions between particles. It this property (like in the SK 
model of spin-glasses) which makes the model solvable in terms of the mean-field approach. Second, 
the particles in this model are point-like. Although formally, due to the proper normalization of 
the parameters of the Hamiltonian, the mean-field thermodynamic limit in this model is quite 
reasonable (see calculations of Section 3), nevertheless, from the physical point of view it looks 
very strange: here we have an infinite number of particles which are confined inside a finite volume. 
Presumably this last drawback could be cured by introducing a finite size for the particles involved. 
On one hand, this should not ruin the mean-field nature of the model, while on the other hand, it 
would make the thermodynamic limit much more reasonable, where both the number of particles 
and the volume of the system are taken infinite while the density of the particles remains finite. 
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